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I The paper by M. Valluri and D. Merritt is discussing an interesting and 

£N| ■ difficult topic, and obviously many questions remain open. I hope, my com- 

ments will be helpful to understand the reasons of the mentioned "debates". 

1. It is not precise to state that "exponential instability of N-body prob- 
lem was shown by Miller [1], Gurzadyan and Savvidy [2]..". The exponents 
in [1] and [2] are different quantities and have different content. 

Miller [1] in his pioneering numerical work has shown the exponential 
growth of separation of coordinates and velocities for N-body systems. It 
was indeed a remarkable work of early period of computer astronomy. Due 
to the development of chaotic dynamics during these decades, our present 
understanding of that problem is much deeper. 

In [2] we have analytically proved the exponential instability of spherical 
N-body systems in Riemannian space with metric depending on the interac- 
tion potential, and the positivity of KS-entropy determined by the curvature 
^ ■ (the curvature was zero in Miller's case). This, according to the theorems 

of ergodic theory, means the existence of the property of mixing and hence 
tending to an equilibrium (local) . KS-entropy is thus defining the time scale 
■ of the tending to equilibrium for large N limit, i.e. the relaxation time; for 

Ci ■ discussion of the meaning of relaxation time see [3]. 

Mixing therefore is not just "a useful way of thinking" but a strictly 
defined property of dynamical systems, with sufficient condition (for typical 
automorphisms) - the positivity of KS-entropy. 

Miller's exponent exists even for N=2, i.e. for two-body problem which 
is integrable, with time scale obviously equal to the crossing/orbital time. 

While exponent in [2] at N=2 turns to infinity and hence - no mixing, no 
relaxation. The systems of 32 particles shown to be exponentially unstable 
by Miller [1], are not unstable via criterion [2]. Indeed, what is the meaning 
of relaxation time for 24 or 32 particles? 

The time scale obtained numerically by Miller, Valluri and Merritt, as 
well as by other authors, must always be of the order of crossing/dynamical 
time, as actually they had obtained, with no direct relation with relaxation 
time scale. Via such numerical experiments hence, for example, nobody 
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had derived (and cannot derive), for example, the two-body relaxation time 
scale, though nobody seems to be alarmed by it. 

The exponent in [2] describes mixing and hence can be associated with 
relaxation, the exponent in [1] - no. Thus, not every exponential instability 
means relaxation (for clear example see [4]). The 'debates' quoted by Vallum 
and Merrit are therefore due to misunderstanding of this crucial point. 

Meanwhile, the relaxation time scale obtained in [2] is supported by: 

(a) observational data on globular clusters (see refs in [5]); 

(b) alternative statistical mechanical methods [6]; 

(c) numerical experiments [7]. 

2. The problem of numerical estimation of Lyapunov numbers for N- 
body gravitating systems is itself a non-trivial one, even forgetting that 
they are limiting in time quantities. In [8] a theorem is proved showing the 
non-equivalence of those Lyapunov numbers and their computer images. 

3. The massive central object is known to increase the instability and 
mixing properties of N-body systems [9]. 
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